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Abstract 

Regarding the adjacency matrices of n- vertex graphs and related 
graph Laplacian, we introduce two famiUes of discrete matrix models 
constructed both with the help of the Erdos-Renyi ensemble of ran- 
dom graphs. Corresponding matrix sums represent the characteristic 
functions of the average number of walks and closed walks over the 
random graph. These sums can be considered as discrete analogs of 
the matrix integrals of random matrix theory. 

We study the diagram structure of the cumulant expansions of log- 
arithms of these matrix sums and analyze the limiting expressions as 
n — > cx) in the cases of constant and vanishing edge probabilities. 

Keywords: Graph Laplacian, Erdos-Renyi random graphs, discrete matrix 
models, cumulant expansion, connected diagrams, Lagrange equation 



1 Introduction 

It counts about thirty years from the time when the map enumeration prob- 
lems of theoretical physics were related with the diagram representation of 
the formal cumulant expansions of the form 

^logE{e«^-} = E ^ -^Cum.m, (LI) 

k>l 

where the average is taken with respect to the gaussian measure over the 
space of hermitian A^-dimensional matrices {H^} and Vn is given by Tr ff^, 
or more generally, by a linear combination of such traces with certain degrees 
q > 3 . This relation between maps, diagrams and random matrices observed 
first by t'Hooft is still a source of numerous results that reveal deep links 
between various branches of mathematics and theoretical physics (see e.g. 
papers [2, 6, 22] for the references and reviews of results). Li these studies, 
the nature of the large- A'^ limit of the right-hand side of (1.1) and existence of 
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its asymptotic expansion in degrees of are the problems of the primary 
importance [3, 7, 10]. 

In present paper we examine the diagram structure of a discrete analog 
of (1.1), where the role of Hjsj is played by the 71 y. 71 adjacency matrix Ayi 
of the Erdos-Renyi random graphs with the edge probability Pn- In this 
setting, Vn of (1.1) is replaced by random variables 

n n 

Xi^)= TrAl = J2{Al)u and y(^^ = {Al),j, q>2 (1.2) 

1=1 *j=i 

and this gives two different discrete matrix models related with the powers 
of adjacency matrix and of the Laplace operator on graphs. Variables X 
and Y represent the numbers of close walks of q steps and walks of q steps 
over the graphs, respectively. 

We study the limiting behavior of the cumulants of Xn^ and Yn'^^ in 
the following three asymptotic regimes, when p„ = 0(1), ^ <C ^ 1) or 
Pn = 0{l/n) as 71 ^ oo. Our main result is that in all of these three regimes, 
the following limits exist (cf. (1.1)) 

-L^Cu7nk{gnYji'i^) ^Fj^''\u;), n ^ oo, (1.3) 

with the appropriate choice of the normalizing factors gn- Limiting ex- 
pressions depend on asymptotic regime indicated by uj G {1,2,3}. The 
same statement as (1.3) holds for random variables Xn^ with corresponding 
changes of gn and the limiting expressions. 

To prove convergence (1.3), we develop a diagram technique similar to 
that commonly adopted in the random matrix theory. We show that the 
diagram structure of the cumulants of (1.3) is closely related with the trees 
with k labeled edges. In the simplest case of q = 2 and co = 2, we derive 
explicit recurrent relations that determine the limits (1.3) and prove that 
the corresponding exponential generating function verifies the Lagrange (or 
Polya) equation. Regarding other values of q and u, we obtain several 
generalizations of this equation. 

Using (1.3), we show that the Central Limit Theorem is valid for the 
centered and renormalized variables Xn^ and Yn'^\ This describes asymp- 
totic properties of the averaged numbers of walks and closed walks over the 
random graphs. Our results imply CLT for the normalized spectral measure 
of the adjacency matrices of Erdos-Renyi random graphs. This improves 
known results about the convergence of the normalized spectral measure 
[1, 16]. Also, we indicate an asymptotic regime when the CLT does not hold 
for variables Xn^ and Yn'^\ 

The paper is organized as follows. In Section 2 we determine two families 
of matrix sums that can be called the Erdos-Renyi matrix models. The 
key observation here is that the graph Laplacian generates the Erdos-Renyi 
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measure on graphs. This allows us to identify as a natural analog of the 
quartic potential V/v = Tr-f^Ar of (1.1). In Section 3, we develop a general 
diagram technique to study the cumulant expansions of the form (1.1) for 
discrete Erdos-Rcnyi matrix models of the first and the second type (1.2). 
In Section 4, we prove convergence of normalized cumulants in three main 
asymptotic regimes. The issue of the Central Limit Theorem is discussed 
at the end of Section 4. In Section 5, we study the classes of connected 
diagrams and derive recurrent relations for their numbers. In Section 6, 
we complete the study of limiting expressions Flf\uj) (1.3) and consider the 
formal limiting transition for the free energy. This free energy takes different 
forms in dependence on the edge probability as n — ^ oo. Section 7 contains 
a summary of our results. 

2 Discrete Erdos-Renyi matrix models 

In (1.1), one can rewrite mathematical expectation as follows 

E jexp i^j^TrH%y^ = ^j^ exp Hj, + Tr cLHm, /? > 

(2.1) 

where iifjv is N x N hermitian matrix, Cjv = C]sf{P) is a normalization 
constant, and the matrix integral runs over the space Hn of all hermitian 
matrices with respect to the Lebesgue measure diJjv- The ensemble of ran- 
dom matrices {H]\f} distributed according to the gaussian measure with the 
density exp{— /3 Tr H^} is known as the Gaussian Unitary invariant En- 
semble abbreviated as GUE . This ensemble plays a fundamental role in the 
random matrix theory (see monograph [19] and references therein). The 
matrix integral of (2.1) is known as the partition function of the matrix 
model with the potential Vn = TrHf^. 

Regarding Ti H% of (2.1) as a formal kinetic energy term, one can 
speculate about the trace of the Laplace operator; being restricted to the 
space of functions on graphs, it leads immediately to a discrete analog of the 
matrix integral of (2.1). In this case the integral of (2.1) is replaced by 
the sum over the set of all n-dimensional adjacency matrices of graphs. This 
approach is fairly natural and one benefits from two important counterparts 
of it. 

From one hand, the use of graph Laplacian indicates a natural analog 
of the matrix models with quartic potentials and clarify relations between 
the weights generated by Tr and TrH^. This helps to distinguish two 
matrix models, the traditional one given by Xn^ and the new one determined 
by Yri'^ (1.2). 

From another hand, a simple but important reasoning relates the graph 
Laplacian with the Erdos-Renyi probability measure on graphs. On this 
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way, the gaussian measure of GUE in the average of the left-hand side of 
(2.1) is replaced by a measure with nice properties we are going to describe. 



2.1 Graph Laplacian and Erdos-Renyi random graphs 

Given a finite graph with the set of n labeled vertices Vn = {^^i, • • • ,Vn} 
and the set of simple non-oriented edges E„i = {e^^\ ■ ■ ■ , e^™)}, the discrete 
analog of the Laplace operator A (7) on the graph 7 is defined by relation 

A(7) = d*d, (2.2) 

where d is the difference operator determined on the space of complex func- 
tions on vertices Vn — ^ C and d* is its conjugate determined on the space 
of complex functions on edges Em C (see [21] for more details). 

It can be easily shown that in the canonical basis, the linear operator 
A (7) = An has n X n matrix with the elements 

r deg{vj), if i = j, 
^ij = S -1' if ^ 7^ j and {vi,Vj) G (2.3) 

L 0, otherwise, 

where deg(t;) is the vertex degree. If one considers the n x n adjacency 
matrix A = A^j) of the graph 7, 

^^fl, if {vi,Vj) e E, i j, 
\ 0, otherwise, 

then one can rewrite the definition of A (2.2) in the form 

n 

Aij = Bij - Aij with Bij = Sij ^ An, (2.4) 

1=1 

where Sij is the Kronecker (5-symbol 



1, if ^ = j, 

0, if iy^j. 



It follows from (2.2) that A(7„) has positive eigenvalues. 

Let us consider the set Tn of all possible simple non-oriented graphs 7^ 
with the set V = Vn of n labeled vertices. Obviously, |r„| = 2"'("'~^)/^. 
Given an element 7 G Fji, it is natural to consider the trace Tr A(7) as the 
"kinetic energy" of the graph 7. Then we can assign to each graph 7„ the 
Gibbs weight exp{— /3 Tr A(7„)}, (3 > and introduce the discrete analog 
of the integral (2.1) by relation 

Zn{P,Q)= E exp{-/3TVAn + Q(7n)}, (2.5) 
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where A„ = A(7„) and Q is an application: r„ ^ R that we specify later. 

Let us note that one should normalize the sum (2.5) by |r„|, but this 
does not play any role with respect to further results. In what follows, we 
omit subscript n in A„. Relation (2.4) implies that 

n n 

TV A = ^ = ^ Ai,- = 2 J2 (2-6) 

i=l *J=1 l<i<j<n 

Then we can rewrite (2.5) in the form 

^n(/3,g)= E e'^^^"^ n e"'''^"- (2-7) 

7ner„ l<i<j<n 

It is easy to see that 

Z„(/3,0) = (l + e-2^)"^""'^^'. (2.8) 

Then the normalized partition function can be represented as 

ZniP,Q) = Z„(/3,Q)/Z„(/3,0) = {e'^W} ^ (2.9) 

where E^{-} denotes the mathematical expectation with respect to the mea- 
sure supported on the set r„. This measure assigns to each element 7 G r„ 
probability 



e-2/3|i?(7)l 



(1 + g-2/3)"("-l)/2 ' 

where £'(7) denotes the set of edges of the graph 7. 

Given a couple i,j G {!,..., n}, one can determine a random 

variable aij on the probability space (r„,P„) that is the indicator function 
of the edge {vi, vj) 

'^^^^ \0, otherwise. 

It is easy to show that the random variables {aij}i<i<j<n are jointly inde- 
pendent and are of the same Bernoulli distribution depending on /? such 
that 

1 - Sij, with probability ^ ^ ^__2$ = (2. 10) 

0, with probability 1 — p. 

The term 1 — 5ij reflects the property that graphs 7 have no loops. 

The probability space (Tji, Pn) is known as the Erdos-Renyi (or Bernoulli) 
ensemble of random graphs with the edge probability p [11]. Since the se- 
ries of pioneering papers by Erdos and Renyi, the as3niiptotic properties of 
graphs (r„,P„), such as the size and the number of connected components, 
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the maximal and minimal vertex degree and many others, are extensively 
studied (see [4, 11]). Spectral properties of corresponding random matrices 
A (2.3) and A (2.4) arc considered in a series of papers (in particular, see 
[1, 13, 14, 15, 16, 20, 24]). In present paper we study the random graph 
ensemble (r„, P„) from another point of view motivated by the asymptotic 
behavior of partition functions (2.9). 

2.2 Quairtic potential and Erdos-Renyi matrix models 

Let us determine the discrete analog of the integral (2.1) with quartic poten- 
tial TiHfj. Once Tr H"^ replaced by Tr{d*d) = TrA (2.2), it is natural 
to consider 

Tr{d*dd*d) = TrA^ 
as the analog of TrH^. Then the partition function (2.5) reads as 

Zn{P,9)= E exp{-/3 TrA„ + 5„ TrA^}, (2.11) 

where gn is to be specified. It follows from (2.3) and (2.4) that 

n n 

Tr A^ = TrB^+ Tr A'^ = ^ {A% + ^ Aj- 

i,j=l i,j=i 

Then, using (2.6) and repeating computations of (2.7) and (2.8), we obtain 
representation 

/l + e-2/3'\"("-i)/2 
ZniP,9n) = ZniP,9n)/Zn{P,0) = f E^,{e^"^"}, (2.12) 

In this relation, we have denoted (3' = (3 — Qn and introduced variable 

n 

^n= auaij, (2.13) 

where Qij are jointly independent random variables of the law (2.10) with 
P replaced by f3'. The average E^/ denotes the corresponding mathematical 
expectation. In what follows, we omit the subscripts P and when they 
are not necessary. 

On can generalize (2.11) and consider mathematical expectation 

n 

E{exp(5„yi''))}, with fW=^(^^),,-, (2.14) 

as the normalized partition function of the discrete Erdos-Renyi matrix 
model that we will refer to as g-step walks model. Also one can consider the 
average 

E{exp(5„4'?))}, with 4«)='IV(^5), (2.15) 
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that we relate with the discrete Erdos-Renyi model for g-step closed walks. 
These models are different analogs of the matrix integral (2.1). 

Finally, let us point out one more analogy between the discrete model 
(2.11) and the gaussian matrix integrals (2.1). We mean the invariance 
property of the probability measure with respect to a group of the space 
transformations [19]. 

Let An denotes the set of of all n-dimensional symmetric matrices whose 
elements are equal to or 1 and the diagonal elements are zeros. It is not 
hard to see that if n x n orthogonal matrix T is such that TAT~^ G An for 
all A G An, then T verifies the following properties 

a) all elements of T take vahics or 1; 

b) any given line of T contains only one non-zero element and 

c) any given column of T contains only one non-zero element. 
Clearly, the set of all such orthogonal matrices is in one-to-one corre- 
spondence with the symmetric group of permutations 5„. 

It is easy to see that T € Xt determines a basis change in because T 
re-enumerates the vectors of the canonical basis associated with the graph. 
Equality Tr (T AT^) = Tr A shows that the probability measure P on An 

P{A) = C-^exp{-/? TrA} 

is invariant with respect the group of transformations 

P(Ty4T^) = P{A). (2.16) 

It is not hard to prove the inverse statement: if the probability measure 
P on the set An verifies (2.16) for all T G 3^„, then the random variables 
given by the matrix elements of A are jointly independent. 

This is in complete analogy with the well-known fact of random matrix 
theory that the invariant gaussian distribution on the space of hermitian (or 
real symmetric matrices) generates independent random variables [19]. We 
do not pursue this topic here and return to asymptotic behavior of the sums 
(2.11) in the limit of infinite n. 

3 Cumulant expansions and connected diagrams 

In this section we develop a diagram technique to study the terms of the 
expansion 

k 

logE {e^^-} = E If (3-1) 

k>i 

where Vn represents represents Xn (2.15) or Yn (2.14). Let us stress that the 
random variables bounded and therefore the series (3.1) is absolutely 

convergent for any finite n and sufficiently small g. 
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Given a vector a = (ii, . . . ,ir), we introduce random variables Ua 

^ ^ r tti^i^Qi^i.^ ■ ■ ■ Qi^i^, r = q for X-modcl 

" \ tti^i^ai^i^ ■ ■ ■ ai^i^+i , r = q + l for y-model 

such that Vn = ^{q}^* ^ai where the sum runs over all possible values of 
is, s = 1, . . . , r. Then one can write relation 

Cumk{Vn) = ^ ■■■ Cum{Uaj_, UaJ, (3.3) 

where 

Cum{Uai,.. ■,UaJ = t^^^ ■ ■ ■ logE{exp(zi;7a, + . . . + ZkUaJ}\zi=0- 

Coefficient Cum{Uai, • • • , f^a^} = Cum{l{{dk)}, ctk = {ai, ■ ■ ■ ,ak) is known 
as the semi-invariant of the family U{ak) = {Uaj}j=i [18]. 

To simplify (3.3), we separate the set {!,... ,n}®('''^) into the classes of 
equivalence according to the properties of the family U(dk)- The rule is 
that given dk, we pay major attention not to the values of variables ij but 
rather to the presence of copies of the same random variable a in U^d^). 
This approach is fairly common in random matrix theory [9, 27]. It leads 
to the diagram representation of the classes of equivalence (see e.g. [6] for 
the review and for mathematical description). This method is used to study 
random matrices with gaussian or centered random variables. We modify it 
to the study of cumulants (3.3) of the Erdos-Renyi models. 

3.1 Connected diagrams 

The diagrams we construct for X-model and y-model are very similar and 
we describe them in common. Let us consider a graph A with r labeled 
vertices {^i,...,^^} with r determined by (3.2). The graph A contains q 
edges: these are ej = {9j,6j-i-i), j = l,...,q; for the X-model we have 
9q-^-l = 9i. Given a = (ii, . . . , ir), one can assign to each 6j the value of ij, 
j = 1, . . . ,r and then each edge Sj denotes a random variable ai-^i-^^. 

To study cum{U{ak)} with given a^, we consider the set Aj. = {Ai, . . . , Ajt} 
of k labeled graphs A; we will say that A; is the element number I of the 
diagram we construct and denote by {e'pYj=i the edges of this element. The 
diagram consists of k elements and a number of arcs a that join some of 
the edges of A's. We draw an arc crj''; that joins ^ and ej', ^ when in 

Here we assume that I < I'. It follows from (3.4) that the arc can have on of 
the two orientations that we call the direct and the inverse one, respectively. 
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li I = I', we assume j < j' and again consider the direct and the inverse 
arcs. 

We say that the edges e^'^ and e^j, ■* represent the feet of the arc aj'i ^ 

and denote this by relation Sj^ ~ Sj', ^ . Clearly, this relation separates the 
set of edges 

£{k,r) = {ef\j = l,...,r;l = l,...,k} 

into classes of equivalence that wc call the color groups of edges. It is 
possible that a class contains one edge only. In this case we say that there 
is a simple color group. Certainly, we color different classes in different 
colors. 




^ K X, ^3 

Figure 1 : A graph A and a diagram 63 for X-model with q = 7. 

An important remark is that it is sufficient to consider the reduced 

diagrams; this means that if there arc three edges that belong to the same 

(I) d') d'') f?' I') 

color group, Ej ~ Ej, ~ £j„ , with I <l' < I", then we draw the arcs a^i ' 

and cr^/;/ •* between the nearest neighbors e^P,e^V and £^!',\e%^ only. The 

same concerns the arcs 0"^^' whose feet belong both to A;. Regarding a 
reduced diagram, each color group has its minimal edge determined in the 
obvious manner. 

To summarize, we consider the set and define the application dk — >■ 5{dk) 
of the set of all to the set of diagrams by drawing all arcs prescribed by 
dk and reducing of set of arcs obtained to the set of arcs between nearest 
neighbors S(afc). As a result, we get a diagram Sj- = 5{dk) = (Afc,E(afc)) 
(see Figure 1). 

We say that a diagram is non-connected when there exist at least 
two disjoint subsets ti and T2 of {1, . . . , n} such that ri U T2 = {1, . . . , n} 
and there is no arc with one foot in A(ri) = {Xj,j G ri} and another in 
A(r2) = G T2}. If there is no such subsets, then we say that the 

diagram is connected. The following statement is a well-known fact from 
the probability theory. 
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Lemma 3.1 . If ak is such that the corresponding diagram 6k = 
is non- connected, then Cum{U{dk)} = 0. 

Proof. It is clear if there exist two subsets ri and T2 as described be- 
fore, then the u-algebras generated by random variables {U^,iJ, G Ai} and 
{Ui,, V G A2} are independent. The fundamental property of semi-invariants 
is that in this case Cnm{W(ajt)} vanishes [18]. 

3.2 Cumulants and sums over connected diagrams 

It is clear that for any given set = {ai,...,afe} there exists only one 
diagram = (5(5^). We agree that two diagrams b^, = (A^j'^k) aiid ^'k = 
(A^, S^) are not equal, 5k / 5^ if in obvious bijection J{Kk) = A^ preserving 
orderings, we have i7(Sfc) 7^ E^. 

Let us say that two vectors q^. and a^. are equivalent, ~ if 
<^(<3fc) = ^{(^'k)- Relation ~ separates the set {1, . . . ,n}®*^'''^^ into the classes 
of equivalence that we denote by C{Sk),6k = S{ak)- Let M{5k) = \C{5k)\ be 
the cardinality of the equivalence class C{Sk)- 

Relation (3.4) means that two edges e and e' belong to one color group 
if and only if random variables a and a' assigned to these edges by are 
equal. This is true for any given G C{6k)- Let us denote by m{5k) the 
number of color groups of 5k- Clearly, random variables that belong to 
different groups are jointly independent. 

Lemma 3.2. The right-hand side of relation (3.3) can he represented as 
■■■Y. Cum{Uo,,,...,Ua,}= Y ^ihWiSk), (3.5) 

where Vk is the set of all possible connected reduced diagrams of the form 
(Ayt, E) and 

W{5k)= (Ea)'"^'''=) (Ea)^(''-^'=) . (3.6) 

In this relation vr^ denotes a partition vTs = (ri, . . . r^) of the set {1, 2, . . . , A;} 
into s subsets; Ilk is the set of all possible partitions, and xi'^s^^k) is the 
number of additional color groups generated by tTs ■ 

Proof. By definition of the semi-invariant [18], we have that 
Cum{Ua,,...,UaJ= ^ (-1)^(5-1)! E{U(Ti)}---E{U(r,)}, (3.7) 

where we denoted U(Tj) = H^Grj ^f^- R easy to see that that the right- 
hand side of (3.7) depends on 6k only, because it does not change when dk 
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is replaced by another a'p, from the same equivalence class C{5k)- Denoting 
by W{5k) the right-hand side of (3.7), we get (3.5). 

To prove (3.6), let us choose an element a^*^^ of the equivalence class 
C((5fc) and consider first the trivial partition ttq = (tq) with tq = {1, . . . , A;}. 
Then xij^si ^k) = ^ and 

where we have use the independence of random variables {dij} that belong 
to different color groups of 5k and the fact that a G {0, 1}. 

Now let us consider a partition vr^ / ttq. It is clear that there exist at 
least one color group of two or more edges such that its elements belong 
to two or more different subsets A(rj) = {A^, t G Tj}. In this case we say 
that this color group is separated into subgroups; each subgroup contains 
edges of this color that belong to one subset Tj. We color the edges of each 
new subgroup in the same color, the edges of the different subgroups - into 
different colors. It follows from (3.7) that the color group separated by tTs 
into v + 1 new subgroups provides the factor (Ea)''"'"^ to the right-hand side 
of (3.6). We say that v is the number of additional color groups generated 
by TTg. Regarding all initial color groups, we get factor (Ea)"*('^'=) (Ea)^" 
with Y^v = x(7rs, 5k)- Lemma 3.2 is proved. 

To study J\f{5k), let us consider the set of vertices of A^, 
@k = {0f\l = l,...,k; j = l,...,r}. 

Given dk, we draw the arcs according to the rule (3.4) and get 5k = 5k{dk)- 
Let us say that equality i^!'^ = i^t ^ of (3.4) identifies corresponding vertices 
9j^^ and 9^', ^ and denote this by 

9f'^9f,\ (3.8) 

It is easy to see that relation (3.8) separates 0^ into classes of equivalence. 
Wc determine the minimal element of the class as the vertex whose numbers 
j and I take minimal values among those of vertices that belong to this 
class. If there is no 9^-, ^ such that 9^'^ = 9j'i \ then we say that 9j'^ belongs 
to the class of equivalence consisting of one element. Let us denote the total 
number of the classes of equivalence by i'{5k)- 

Lemma 3.3. Given 5k G Vk the cardinality Af{5k) = \C{5k)\ verifies 
asymptotic relation 

Ar(^,) = n^(5fc)(i + o(l)), (3.9) 

in the limit n — ^ oo . 
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Proof. It is clear that two equivalent vectors ak ~ generate the same 
partition of into classes of equivalence. Inversely, given a diagram and 
regarding corresponding partition of Qk, we get a separation of the set of 
variables 

Ik = {i'j\ j = l,...,r; l = l,...,k} 

into viSk) groups. Variables that belong to the same group are equal between 
them. To get all possible dk from the same equivalence class, we allow 
variables to take all possible values from 1 to n with obvious restriction 
that variables from different groups take different values. Then obviously 

J\f{Sk) = n{n — as n — oo. (3.10) 

Lemma 3.3 is proved. 

Let us complete this subsection with the following useful remark. The 
set of diagrams we constructed gives a graphical representation of vec- 
tors dk and describes the classes of equivalence of such vectors. One can 
push forward this representation and consider the set of graphs of diagrams 
Gk = G{5k) generated in natural way by diagrams 5k by gluing the edges of 
elements A in the way prescribed by the arcs of 5k- Actually, this is what is 
usually done in the standard diagram approach of random matrix theory. 

The graph G{5k) = (B,^) has i'{5k) labeled vertices that correspond 
to the minimal elements of the classes of equivalence of Gj^. The vertices 9 
and 9' are joined by an edge 9') G ^ if there is an edge e G £{k, r) that 
joins corresponding classes of equivalence of Qk- Certainly, the number of 
edges of G{5k) is equal to the number of all color groups m{5k)- 

The graph representation is very useful when 5k contains A; — 1 arcs 
only. We say that 5k are tree-like because in this case, the graphs G{5k) of 
y-model are trees with color edges. 




Figure 2: Tree-like diagram and corresponding tree for y-model with g = 3 

Remembering that the arcs cr G S have the direct and the inverse ori- 
entation (3.4), we agree to consider the graphs G(6k) with the edges of 5k 
glued in the inverse sense only. In what follow, we will use representations 
of dk by diagrams 5k and graphs G{5k) both. Note that '^{5k) has no loops. 
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4 Diagrams and limits of the cumulants 



In this section wc characterize the classes of connected diagrams that provide 
the leading contribution to the cumulants of X and Y models. The following 
three asymptotic regimes known from the spectral theory of random matrices 
are also distinguished in present context. We refer to them as to 

• the "full random graphs" regime, when p = Const as n — oo; 

• the "'dilute" regime, when p = Cn/n and 1 ^ c„ <C n; 

• the "sparse" regime, when p = c/n and c = const as n — > oo. 

One can consider also the fourth asymptotic regime of the "very sparse" 
random graphs, when p = c^jn and c„ — > as n — > oo. In this regime, 
the diagrams of the leading contribution to cumulants are degenerated, so 
we formulate corresponding results as remarks to the proofs of the main 
theorems. We omit subscripts in Cn when they are not necessary. 

4.1 Asymptotic behavior of cumulants 

In this subsection we formulate the main results on the asymptotic behavior 
of the cumulants Cum^iTJ^. We start with the Y-model. In what follows, 
we denote by Un x Vn asymptotic relation n„ = + o(l)) as n — ^ oo. 

Theorem 4.1 Asymptotic behavior of Cumk{Yn''^) is determined by the 
following classes of diagrams: 

A) in the full random graphs regime, the leading contribution to the sum 
(3.5) is given by the diagrams V^i^\y) that have exactly k—1 arcs; the graph 
G{6k),6k G 'D^\y) is a tree; then 



B) in the dilute random graphs regime, the leading contribution to (3.5) 



C) in the sparse random graphs regime, the leading contribution to (3.5) 
IS given by the diagrams ^^c/i that the graph G{Sk),Sk € 

is a tree Ti with I edges, 1 < I < {q — l)k + 1; then 

(g-l)fc+l 




(4.1) 







(4.3) 



1=1 5keV(3)(Y),GiSk)=Ti 
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Corollary of Theorem 4-1- Remembering that Ea = p = Pn, we can 
reformulate the results of Theorem 4.1 in the following form: 

lim -1^ Currik ( , Y^A = F^'^^ (co), a; = 1,2, 3, (4.4) 

where the numbers pjf^ (i) represent the contributions of corresponding fam- 
ilies of diagrams (weighted by W{Sk) in the full random graphs regime). We 
discuss relation (4.4) in more details in subsection 4.3. Explicit expressions 
for some of F^\u) will be obtained in Section 6. 

Let us consider the X-model. Now the difference between the cases of 
the even and odd numbers q becomes crucial. 

Theorem 4.2 Asymptotic behavior of Cumk{Xn) is determined by the 
following classes of diagrams: 

A) in the full random, graphs regime, the leading contribution to the sum 
(3.5) is given by the diagrams 'D^\x^'i^) that have exactly k — 1 arcs; in 
this case 

C7nmfc(x('?))x Yl W^(4)nM'=+^ (4.5) 

where expressions for W(6k) is given by (3.6) with 6k G 'D^j^\Xn ''); 

B) in the dilute random graphs regime, the leading contribution to (3.5) 
in the case of X -model with even q = 2q'is given by the diagrams (X) 
such that the graph G{5',^),d'f. £ (X) is a tree with the maximal possible 
number of edges; then 

Cumfc(42.')) X (Ea)('''-i)*^+i n^'i' -^)^+-' ■ (4.6) 

5,ei5f {x(V)) 

in the case of X -model with odd q the leading contribution to (3.5) is given 
by the diagrams {X) such that the graph G{5'j^),5'j^ G {X) is a cycle 
with he maximal possible number of edges; then 

Cumfc(X(2?'+i)) X J2 (Ea)2<?'+1 n^^'+i; (4.7) 

C) in the sparse random graphs regime, the leading contribution to (3.5) 
in the case of even q = 2q' + 1 is given by the diagrams V^^^ (X) such that 
the graph G{6'^),d'^ G V^^\X) is a tree Ti with I edges, 1 < I < (q' - 1)A; + 1; 
then 

(q'-l)fc+l 

CurukiXi^'i'^) ^ J2 E (Ea)'n'+^ (4.8) 

1=1 5fc6©(3)(X),G(5fc)=T; 
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in the case of X-model with odd q = 2q' + 1, the leading contribution to 
(3.5) is given by the diagrams T)^^\X) such that G{6k) is a graph T^^^ with 
I edges and one cycle only; then 

(g'-l)fe+l 

Cumfc(42.'+i)) X Yl E (4-9) 



Corollary of Theorem 4- 2- In the full random graphs regime, 

lim ^ Cumk ( r ^^i^^) = (4-10) 

and in the dilute and sparse regimes we have, respectively, 

Cumu {^,X^^''^) = = 2,3, (4.11) 

and 



lim ^ 

Kc<n cn \ci 



and 



,4!S„ ^Cum,(X<V+i)) ^ a>W+l)(2, (4.12) 

We discuss relations (4.10)-(4.13) in Subsection 4.3. Explicit expressions for 
some of (^^^\(jo) will be obtained in Section 6. 

4.2 Proof of Theorems 4.1 and 4.2 

Proof of Theorem 4-1- In the full random graphs regime, Eo is a constant 
and all the terms of the sum (3.6) are of the same order of magnitude. The 
leading contribution to (3.5) comes from the diagrams 5k such that the graph 

G{5k) has maximally possible number of vertices. In this case the number of 
arcs in 5k is minimal, i.e. is equal to A: — 1 and G{5k) is a tree. This proves 
(4.1). 

In the case of dilute and sparse graphs, relation Ea^ = o(Ea) shows that 
the leading contribution to the sum (3.6) is obtained from those partitions 
TTs that x{'^s'i5k) = 0. Since 5k is connected, then only trivial partition ttq 
verifies this condition. Let us denote by I the number of edges in G{5k)- If 
G{5k) is given by atree, then v{5k) = 1 + 1 and {Eaf\M{5k)\ = nc'(l + o(l)). 
If G{5k) is not a tree, then by the Euler theorem, v{5k) < I and such diagrams 
provide a contribution to (3.5) of the order o{nd'). In the case of dilute 
random graphs, c — oo and the leading contribution is obtained from the 
graphs withZ = [q — Vjk + l; other trees do not contribute. This proves (4.2). 
To show (4.3), it remains to note that in the sparse random graphs regime 
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all trees with 1 < I < (q — l)k + 1 provide contributions of the same order 
of magnitude. Theorem 4.1 is proved. 

Proof of Theorem 4-2- In the full random graph regime, the leading 
contribution to (3.5) comes from those 5^ that have fc — 1 arcs. This implies 
relation (4.5). 

It is easy to see that in the dilute and sparse regimes, the only trivial 
partition ttq contributes to (3.6). In the case we repeat arguments 

of the proof of Theorem 4.1 and get relations (4.6) and (4.8). 

Let us pass to the case of odd q = 2q' + 1 and consider diagrams Si with 
one clement Ai. It is clear that graphs G{Si) always contain at least one 
cycle. Then G{6k) also contain at least one cycle. If the number of edges 
of G{5k) is equal to then the contribution of such a diagram is of the 
order (Ea)'n'+^~^, where w is the number of cycles in G{5}.). The graphs 
with one cycle only provide the leading contribution of the order (Ea)'n'. 
In the dilute random graphs regime the leading contribution comes from 
the graphs with maximal number of edges. Then we conclude that in this 
case G{8k) represents a cycle with 2q' + 1 edges and 2q' + 1 vertices. This 
proves (4.7). In the case of sparse regime, we consider 5^ such that G{5k) 
has exactly one cycle. This implies (4.9). Theorem 4.2 is proved. 

Remark. It follows from the proof of theorems 4.1 and 4.2 that if Pn = 
c/n and c ^ as n — oo, then the leading contribution to Cumfe(Y^^^) and 
Cumk{Yn'^) is given by the connected diagrams 5k such that G{Sk) have 
minimal number of edges. It is easy to see that in this case 



Cumk{Y(i^) = 2^-^pnn\l + o(l)) = 2^-^cn{\ + o(l)), (4.14) 



CMmfc(X^2y')) = 2^-^pnn^{l + o(l)) = 2''-^cn{l + o(l)). (4.15) 
Also we have Cumk{x'i'^'^^^) = 0{pln^) = 0{c^) as n ^ oo, c ^ 0. 



4.3 Central Limit Theorem for variables X and Y 

Regarding the definition of the cumulants, it is easy to see that for any 
constant C, 



and 



Cumk{Vn + C) = CurukiVn) + C6k,i. 
Also, if there exists such a sequence 6„ that 

^CunikiVn) (pk, k> 1, 



(4.16) 



'n 



then 




(4.17) 
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Relation (4.17) means that the probabihty law of the centered random vari- 
able Vn = {Vn — EV^)/-v/5^ converges to the Gaussian distribution A/'(0, 02). 
Regarding centered and normalized variables Xn^ and Yn'^\ we formulate 
the Central Limit Theorem. 

Theorem 4.3 If n ^ oo, and Pn determines one of the three asymptotic 
regimes indicated by lo, the following random variables converge in law to 
the Gaussian distribution: 

" " -^m4'\^)), = 1,2,3, (4.18) 



Vp^ (p„n)9-i 



and 



(2?') _ 17 v(29') 

M{0,^^2'''\io)), w = 2,3. (4.19) 



Also convergence to the standard normal distribution holds in the case of X 
model in the full random graphs asymptotic regime; 

— ^ ^ ^ ^^^^ ^(,)^^ ^^ 20) 



and in the regime of dilute random graphs for X^'^^ with odd q = 2q' + 1: 

- EX^«'+i)) ^ AA(0, $?'^'+^)(S)). (4.21) 

There is no convergence to the normal distribution of Xn'^ in the regime 
of sparse random graphs. 

Proof. The proof follows immediately from relations (4.4) and (4.11) 
with properties (4.16) and (4.17) taken into account. 

Let us discuss relations of Theorem 4.3 with the spectral theory of ran- 
dom matrices. Remembering that PnU = c and introducing the normalized 
adjacency matrices 

^(n,c) ~ "^y^-An, 

we can introduce variables Mi",'*^-* = - Ti A'^'^ ^ and rewrite (4.19) in the 

2q' n (n,c) ^ 

form 

V^c{M^y^-BM^y^)^M{0,^f'''\u)), a; = 2,3. (4.22) 
Also, we deduce from relations (4.11)-(4.13) that 

}}^™2T^ = ^?'''H^) = m2A^) and JimEM§'^\ = 0. (4.23) 
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Variables Mq^''^^ represent the moments of the normahzed eigenvalue 
counting measure of random matrices A(„ (,)! 

^n(A) = #{j : Xf'"^ < X}n-\ 

Convergence (4.23) implies the weak convergence of the measures an- This 
convergence is established in the present and more general settings in papers 
[1. 13] and [16]. Central Limit Theorem (4.22) improves results of these 
papers. 

Let us note that relations (4.14) and (4.15) imply that 

^CumkiY^'i^) ^ 2^-^ as n^oo,c^O. (4.24) 

The same is true for the cumulants of Xn'^ ■* . This means that the Central 
Limit Theorem does not hold for these variables in the asymptotic regime 
of very sparse random graphs. 

Returning to variables Yn*^^ , we rewrite (4.4) with A; = 1 in the form 

lim — Ey(9) = Fi'^\u) = 1, (4.25) 

since there exists only one connected diagram. Variable Yn'^^ counts the 
number of g-step walks over the graph, and relation (4.25) shows that this 
number is asymptotically proportional to nc'^. It is natural to expect this 
result because the average degree of a vertex in the Erdos-Renyi random 
graph converges to c, as n — > oo [11]. It follows from (4.18) that 



7n 



nci " \nci ^ J 

where 7„ converges in law to a gaussian random variable. Then convergence 
with probability 1 holds 

— Fi-?) ^ 1 
ncfl " 

as n ^ oo in the dilute and sparse random graphs regimes. The limiting 
variance of 7„ depends on the asymptotic regime. We return to this question 
in Section 6. 



5 Number of tree diagrams 

In this section we derive recurrent relations that determine the number of 
connected diagrams 5^ = (A^, S) with minimal number of arcs, 



p(') = {(Afe,S): |S| = /c-l}. 
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In the case of y-model, the graphs G{6k) have the tree structure and we 
refer to Sk as to the tree diagrams. In the case of X model we refer to such 
diagrams as to the tree-like diagrams. It is clear that 

|p(^)(x('?))| = |p«(y('?))|=4^\ 

so we consider the case of the F-model only. 



5.1 Recurrent relations and Lagrange equation 

First let us remind that the arcs in 6^ are interpreted by G{6k) as gluing 
between edges that make a color group. Let us say that the edges that are 
not glued are not colored and stay grey. So, we forget the color of the simple 
color groups consisting of one edge only. The different color groups and grey 
edges correspond to the edges of G{Sk)- It is convenient to color the arcs 
that join the edges of one group in the same color as the edges of this group. 
We start with the following simple statement. 

Lemma 5.1. Let a diagram 6^ G 'D^j^\y'^'^'>) have ares of s different 
colors. Then there are {q — l)k — s + 1 grey edges in 6k- 

Proof. Let us assume that there are jij arcs of each color. Obviously, 
fii + . . . + jjLg = k — 1 and the total number of colored edges is A; — 1 + s. 
Taking into account that the total number of edges in Sk is qk, we obtain 
the result. Lemma is proved. 



Lemma 5.2 Given q > 2, the numbers d), 
are given by equalities 



(<?) 



\vI^\y(i^)\ with k > 1 



d 



(?) 



A;! 



{q-l)k+l 



h 



(q) 

k ' 



(5.1) 



where the sequence {h^k ^}k>o is determined by the recurrent relation 



(,) _ (g-l)fc + l 
~ k 



E 



jl+---+jq=k-l 

ji>0 



\ 



■hf 



, k>l 



(5.2) 



with initial condition h, 



(g) 



1. 



Proof. Let us remind that since our diagrams are reduced, then each 
color group has a unique maximal edge. Now let us count the number of the 
diagrams 6k such that the last element contains one and only one edge 
that is the maximal edge of a color group. This means that only one arc 
ends at by its right foot (or leg). Let us denote this arc by a. 
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Figure 3: Example of 5k with possible positions for for the left foot of cr 

Obviously, one can choose one of the q edges of Afc to put this right leg. 
The left leg of a joins with connected tree diagram 6k-i G 'D^}_^{Y^'^^). 
By Lemma 5.1, there are {q — 1)(A; — 1) — s + 1 grey edges and s maximal 
edges of s color groups, where one can put the left foot of a. Then we obtain 
relation 

Now let us consider the case when there are / arcs ai, . . . ,ai that end 
at Ajfc with 2 <l < min{q,k — 1). There is (f) possibilities to choose the 
emplacements for the right legs of these arcs. To put the left legs of these 
arcs, one has to choose the / subsets A^^^ . . . , A*^') of ji = |A(*)| elements such 
that ji + . . = k — 1 and ji > 1, then to create I connected sub-diagrams 
(A(*), i = and to choose the emplacements for the left legs of 

ai in (A(^),5(*)). This produces 

...+ji=fe-i -^i- j^- i=i 
ji>i 

diagrams. Summing up, we derive for numbers {d^^\k > 1} with given 
q>2 the following recurrent relation 

= E%,^-ii(Ox (^) E ri((.-i).-.+i)4f' (5-3) 

1=1 VV n+...+ji=k-iJ^- 1=1 

ji>l 

with initial condition d!f^ = 1. Here we denoted by /[i,a;-i](0 the indicator 
function of the interval [1, A; — 1]. 

Introducing the auxiliary numbers hj , such that 
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we reduce (5.3) to relation 
(9-l^)fc + l ^^ = E hhk-i] (0 X (^^ Yl ^hhh ■■■K,k>\. (5.4) 

If we introduce an auxiliary number /iq = 1; then we can rewrite (5.4) in 
more compact form. Indeed, we can make an agreement that each labeled 

(k) 

edge ei of the element serves as the right foot of the corresponding 
arc ai but some of these arcs can have the right leg "empty", i.e. with no 
elements attached to their right feet. Then the variable I of the sum of (5.4) 
represents the number of non-empty arcs, and (^) stands for the choice of 
these non-empty arcs from the set cji , . . . , (jq . This corresponds to the choice 
of the variables ji in the product hj^ ■ ■ ■ hj^ that take zero value. Finally, 

we set hj^^ = hj,j > and conclude that relation (5.4) is equivalent to 
recurrent relation (5.3). Lemma 3.5 is proved. 

Lemma 5.3. 

The generating function 

k=0 

is regular in the domain {z : \z\ < l/(g^e)} and verifies there the analog of 
the equation 

H,{z)=cxp{qz{Hg{z)y-'}. (5.6) 

It follows from (5.3) that the numbers h^j^^ can be found from recurrent 
relations 

Jl+---+iq-l=fe 

ji>0 

Remark. Relation (5.6) with = 2 is known as the Lagrange (or Polya) 

equation [17, 23, 25]. In the next subsection we give another derivation of 
(5.6) by using the notion of color trees. In Section 6 we obtain one more 
generalization of the Lagrange equation. 

Proof of Lemma 5.3. It is easy to deduce from (5.2) that Hq(z) verifies 
integral equation 

H,{z) -l = ziq- l){H,{z)y + j\H,{Oy 

Jo 

that is equivalent to the differential equation 

^ = -^(^-(-^^^ H,iO) = l. (5.8) 

dz l-{q^-q)z{H,{z)r^' ' ^ ' 



21 



Substitution 

^{z) = z{H,{z)f-^ (5.9) 
transforms (5.8) to an elementary equation 

1 V(^) 



z 1- (g2 -g)^(z)' 



Resolving this equation with obvious initial condition, we conclude that tpi^z) 
verifies the algebraic equation 

ij{z) = ze^^i-^^^^'^ (5.10) 

known as the Polya equation. Then (5.6) easily follows from (5.9) and (5.10). 

Using the standard method of the contour integration, we deduce from 
(5.10) explicit expressions for the coefficients ipk of the generating func- 
tion ip{z) = Yl,k>o'^kZ^ ■ The first observation is that the inverse function 
ip*{w) = ?x;e~'^('?~^)^ is regular in the vicinity of the origin. By the Cauchy 
formula, we have 

By changing variables hy z = 'ip*{w), we get dz = {l — q{q — l)w)e~'^^'^~^^^dw 
and find that 

27r i / w'^ 



Then 

il,,. = (n(n- 

{k-iy. 



L.fe-2 

^, = (q(q-l)f-^——, (5.11) 



and (5.7) follows. Lemma 5.3 is proved. 
5.2 Color trees 

In this subsection we give another derivation of the Lagrange equation (5.6) 
based on the graph representation G{5k) of the diagrams S^- For simplicity, 
we consider the case of Y-model with q = 2 only. 

If 6k has k — 1 arcs, then G{dh) is a tree. We say that the arcs that 
indicate the groups of edges to be glued are all of the same color. Then the 
corresponding edges of G{5k) are of this color. We say that the edges that 
are not glued are grey. 

Let us derive expressions for the number of the elements of the set 

Tk = IGkl, Qk = {G{6k),6k e P«(y(2))}. 

We denote by the number of rooted trees of the form G{5m)- 

Let us consider left element of Ai as the root edge p for the trees we 
construct (see Figure 4) . We can attach to this root r elements vi,...i'r 
that we choose by {^) ways. We glue these elements by their right edges. 
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(1) 



Figure 4: Rooted color tree and color tree 



Then we separate m — r elements into tlic groups oili,l2-, ■ ■ ■ J, elements 
and construct rooted trees and attach them by their roots to the free ele- 
ments of 1^1 , . . . , . There are possibilities to separate m — r elements 
into these groups. Then we get 

r=l V J l^+...+l^=m-r 

where the sum runs over all > and we accepted that Tq = 1. 

Simplifying (5.12) and taking into account that Tq = 1, we derive from 
(5.12) the following relation for the generating function t{z) = YlikXi^k^'^ i 
tk = fk/k\ (cf (5.6)): 

t{z) = ex.^{2zt{z)}. (5.13) 

Then Tm = 2"^(m + and the numbers tm are given by recurrent 

relations 

Ul \ m— 1 

tm = X! tjtm-l-j, to = 1. (5.14) 

TTl 

j=0 

Now let us return to the construction of the tree G{6k)- Now the left 
edge of Ai can serve as the root for the left rooted tree of mi elements and 
the right element serve as the root for the rooted tree of m2 = A; — 1 — mi 
elements. Then 

^ (fc-1)! ^ ^ 

-'-k — / ^ I I -'■mi -'■7712; 



mi>0 



mi\m2'- 



where Tq = 1. Then 



Tk _ 2*=(fe + l)*^-2 

tjtk-l~j — , j_ -.tk 
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Then equality Tk = 2^{k + if'"^ follows. 



In the general case q > 2, one can easily repeat this procedure of tree 
construction and derive (5.6) directly, without use of the numbers h^^\ 

6 Limits of the cumulants 

In this section we specify the limits of the cumulants of the F-model F^\i) 

(4.4) and of the X-model ^^^-'(i) (4.10)-(4.13). Wc start with the case of 
full random graphs regime, but the most attention is paid to the dilute and 
sparse random graphs regimes. 

6.1 Cumulants of full X and Y models 

In Section 4 we have shown that the leading contribution to the cumulants 
of (4.1) and (4.5) is given by those diagrams 8k that have exactly k — 1 arcs 
and the graph G{5k) with the maximal number of vertices. The weight of 
the diagram W{5k) is determined by relation (3.6), where all terms {Eaf 
are of the same order of magnitude (see the proof of Lemma 3.2). Regarding 
y-model and taking into account that in this case m{5k) = {q~ l)k + l (see 
Lemma 5.1), we combine (3.6) with (4.1) and conclude that (4.4) is true 
with 

Fi'\^)= E E (6-1) 

where x{T^s'-,^k) is determined in Lemma 3.2. 

Considering X-model, it is easy to see that (3.6) gives the same expres- 
sion for W{5k) as for the y-model. Using (4.5), we arrive at the conclusion 
that (4.10) is true with ^^^\A) = fI^\a) (6.1). 

Restricting ourself to the first two cumulants, it is easy to compute that 

F^«)(l) = $(9)(1) = 1 

and that 

6.2 Cumulants of dilute F-model 

Regarding (4.2) with Ea = c/n and using results of Section 5, we conclude 
that (4.4) takes the form of 

^^'^(2) = Jm^ ^Curuk {Jz^y^i'^) = 2^-^|2^?^(r)| = 2^-'d!i\ (6.2) 
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where c?^'^'* are determined by relations (5.1) and (5.2). In (6.2), we have 
taken into account that the A; — 1 axe of the diagrams 5^ can be drawn in 
the direct and inverse sense. This produces the factor 2^~^. 

Let us consider in more details the case of g = 2 that corresponds to 
the continuous matrix model (2.1) with the quartic potential. In this case 
relation (5.2) takes the form 

A2) _k + l^ V .(2), (2) ,(2)_, 

3l+32=k\,ji>0 

It follows from (5.7) that in this case hf^ =2^{k + lf-'^/k\ and 

=2\k + lf-''. (6.3) 

The combinatorial meaning of ' can be explained with the help of the 

notion of the dual diagrams that wc briefly describe below. 

Regarding 5^ £ T)^^'' iY^"^^), we construct the dual diagram dk by trans- 
forming the edges of elements \j into vertices joined by an edge; the arcs 
remain without changes and indicate the vertices of 5k to be glued (see 
Figure 5). It is easy to see that the graph G{5k) is a tree. Then d\. is equal 
to the number of all non-rooted trees constructed with the help of k labeled 
edges. Since the edges of 5k are labeled to be the left and the right one, we 
can think about the edges of G{5k) are the oriented one. 

Regarding the dual diagrams in the general case of g > 2, we see that 
these are trees when q = 2 and trees constructed from oriented chains of q 
edges when q > 2. It is not clear, whether it is possible to derive from (5.7) 
explicit expressions for d^^^ with q > 2. 

6.3 Cumulants of sparse y-model 

The diagrams of V^^^ (Y) are of more complicated structure than those of 
V^^^Y) and we can not obtain recurrent relations to determine Flf\3) in 
the general case of g > 2 and A; > 1. We present the results for Fj^^'*(3) and 
Fi'^\3) only. 
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(2) 

Lemma 6.1 The coefficients (3),/c > 1 are determined by relations 

2^-^{k - 1) 



Ff (3) 



where Wk are determined by relation 



-Wk, 



k—s 



^WjWk-s-j, k>l, 



(6.4) 



(6.5) 



and the nvmbers Wj,j > are such that the generating function W{z) = 
J2k>o WkZ^ verifies equation 



W{z) = 1 - ce^^ + cexp \2z[e'^^W{z) - 1]} . 



(6.6) 



Proof. Following the lines of subsection 5.2, we consider first the set of 
the rooted trees of the type we are interested in and denote their contribution 
by T^, where m is the number of elements A. 





ei 



(1) 



Figure 6: Elements to construct a rooted tree and example of the tree 

In the present case the trees are constructed with the help of elements A, 
but there can be double elements that are glued to one of the edge, and there 
can be simple elements, and there are elements glued to simple elements all 
the long (see Figure 6). Let us consider the root edge p. Denoting by u 
the number of double elements attached to p and by r the number of simple 
elements attached to p by one edge. Also we denote hy Vi, i = 1, . . . , r the 
numbers of elements glued all the long to these r simple elements. Taking 
into account possibilities to choose u, r, and Vi elements from the set of m 
labeled elements, we obtain the recurrent formula 




E 



2^1 . . . 2""^ 



Vl+...+Vr=V,V=l 



m — u 



Vl,. 



26 



, , ]Ti,---Ti^, (6.7) 

h+...+lr=m^u^r-V \ n'---;'?- / 

where we denoted the multinomial coefficients 



m — u — r\ [m — u — r)\ 



\vi, . . . ,Vj. j vi \ ■ ■ ■ \ {m — u — r — V)\ 

Factor c in the right-hand side of (6.7) takes into account the contribution 
(or weight) of the root p. Simplifying (6.7) and denoting = T^/m}. for 
m > 1, we get relation 

"1 nu rri—u ryr m—u—r r\vi nvr 

A V — ^ ^ X — ^ \ — ^ ' ' ' ^ — ^ ^ ^ 



u=0 r-=l ■ y=i;i+...+i;^,\/=l 1" ^" h+...+lr=m-u-r-V 



Passing to the generating function W{z), we derive equality 

W{z) - 1 = ce^^ [exp \2ze^^W{z)^ - l] , (6.8) 

that is equivalent to (6.6). Let us note here that (6.8) generalizes the Polya 
equation (5.13). 

Returning to the set T)j^ \Y), we consider Ai as the simple root element. 
Assuming that there are s — 1 elements glued all the long this root element, 
we dispose of A; — s elements to construct the rooted trees using the left 
and the right edge of Ai as the roots pi and p2- Denoting the number 
Tfe = \Df\Y)\, we get relation 



— 1/ ^ \ mi, 1712 1 

S=l \ / m-i+m2=k-s \ ^' ^/ 

mj>0 



Simplifying this relation and passing to variables Wm = Tm/{ni — 1)!, we 
obtain equality (6.5). This completes the proof of Lemma 6.1. 



In the general case of g > 2 we determine the first cumulant only, 

n 

(3) = lim E E {Ai,i2Ai2i, ■ ■ ■ Ai^i,^, ) . (6.9) 

c=Const «iv,«9+i=l 

The vertices of A are ordered and we denote 6i= p and 62 = v. It follows 
from Theorem 4.1 that f['^\^) is determined by the number of diagrams 5i 

consisting of one element A such that corresponding graph G{Sk) is a tree 
of / vertices, 2 < I < q. This problem is closely related with the studies 
done in [1, 13] and formalized in [16]. These studies are related with the X 
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model and we present the results in the next subsection. In what follows, we 
describe briefly corresponding approach and modify it to suit our situation. 

To study the right-hand side of (6.8), it is useful to introduce the notion 
of a walk ^. A walk (^g of q steps is an ordered sequence of q letters, starting 
with p and followed by v always. Given a sequence /g+i = {ii,i2, ■ ■ ■ , iq+i), 
we construct corresponding = ^(/g+i) by the following recurrent rule: 
regarding the value ig+i, s > 1, wc compare it with the values of {ii, . . . , 1,} 
and write a letter number s + 1 in ^. If there is no < j < s such that 
is+i = ij, the we write a new letter that is not present in ^g- If there is such 
j, 1 < j < s such that Ig+i = ij, then we write on the s + 1-th place the 
letter number j of ^s- 

We can represent the walks graphically, by drawing the root vertex p, 
the next vertex v and corresponding edge and then passing along and 
creating the new vertices at the instants, when the new letters occur in ^g. 
We also draw the new edges in this case. 




6162 p p 

Figure 7: A diagram 5i, the walk and corresponding graph G{Si) = G{^q) 

It is clear that there is one-to-one correspondence between the diagrams 
5i and walks ^g. Also, the walk ^g generates in natural way the graph G{^g) 
isomorphic to G{Sk)- 

Lemma 6.2. The number Fi^\3) = FgC^~'^ is given by the sum 

Fg = CFg_^+YlW^ (6-10) 

r=2 

where the numbers Fg{r),r > 2 are determined by the following recurrent 
relations 

v=2 u=0 s=0 \ I 2 J / V2i ^) 

In (6.11), [x] denotes the largest integer less or equal to x and the following 
initial conditions are assumed: Fiifi) = 61^0, Fr-j{r) = for j > and 
F2{2) = 1. 
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Remark. Relations (6.10) and (6.11) are similar to those derived in [1, 13] 
for the number of walks of even number of steps such that their graph is a 
tree. 

Proof. Let us consider the set Hq(r) of walks such that there are r steps 
that start or end in p. We denote by fq{r) = |Sq(r)| the cardinality of this 
set. 

[ q-u-v 

Figure 8: The first-passage decomposition and example of tlie walk 

Let us consider the subset Sq(r, v, s) C Hg(r) consisting of the walks that 
pass the edge ei = (p, v) v times and that have s steps attached to u that 
do not pass ei (see Figure 7). We denote the set of such steps u ^ t: with 
TT / p by Sy^p. Certainly, such a walk has r — v steps that belong to the 
set Sp^ii. 

Any walk of Hg(r, w, s) can be separated into three parts or in other words 
into the following three subwalks: the first one contains the steps p ^ u and 
u ^ p only, the second one is attached to u with the elements of S^^p (the 
left part of the walk) and the third one is attached to p by the elements of 
Sp^p (the right part of the walk). We denote the left and the right parts 
by Sq_^_y{s) and by S„(r — v), respectively, where u denotes the number 
of steps in the right-hand part of the walk. We denote the by Eq{r,v,s, u) 
the set of elements of Sg(r, v, s) such that the right-hand part of the walk 
consists of u steps. 

Then we can write equality 

\^q{r,v,s,u)\ = (^"yit/^^) I I) fq-u-v{s)fu{r-v), (6.12) 

where we have taken into account obvious equalities |S^_^_^(s)| = |Sq_„_„(s)| 
and |H(^(r — v)\ = |H.u(r — v)\. 

Let us explain the from of the combinatorial coefficients in the right- 
hand side of (6.12). If r and v are both even, then there are v/2 steps in 
St^^p and r/2 — v/2 steps in Sp-.^^. Therefore there is a choice to perform 
each of the r/2 — v/2 steps after one of v/2 steps u ^ p. This produces the 
second binomial coefficient of (6.12). 

If r is odd and v is even, then the last step of the form p ^ it with it ^ p 
is to be performed after that all the steps u ^ p are done. Then we get 
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again the combinatorial coefficient as that given by (6.12). If v is odd, then 
there are stih [v/2] steps u ^ p. 

Let us explain the first binomial coefficient of (6.12). If s and v are both 
even, then there arc C^'^^2-]~^) possibilities to perform s/2 steps of S^^p 
after v/2 steps of Sp^y. If s is odd, then the last step is performed after 
that all v/2 steps are performed. If v is odd, then there are [u/2] + 1 steps 
p^v. 

Summing (6.12) with respect to parameters a,?x, and s, and taking into 
account that Eo, = c/n, wc get relation (6.11) for the contribution to Fq 
determined by the walks with r, r > 2 steps attached to p. If r = 1, then 
corresponding contribution is given by cFq_i. Summing all contributions, 
we obtain formula (6.10). Lemma is proved. 

Remark. It is easy to see from (6.12) that fz{2) = 0. Moreover, it is not 
hard to show that /2/+i(2s) = for all l,s > 1. This is in agreement with 
the obvious observation that the set of walks H2;+i(2, s) is empty. 



6.4 Cumulants of dilute and sparse X-model 

Analysis of cumulants of X-model in the dilute random graphs regime is 
very similar to that of the Y-model in the same regime. 

Lemma 6.4 Relations (4-11) o-nd (4-1^) are true with 



^^^'i'+^\B) = {A(^ + 2f-^ and ^^^'^'\b) 



1 / 2q'^ 



2' +l\q' 



(6.13) 

where d^ is determined by relations (5.1) and (5.2) with q replaced by q' . 

Proof. According to Theorem 4.2, the leading contribution to (4.7) in 

the case of odd q = 2q' + 1 is given by the diagrams 5k such that the graph 
G{6k) is of maximal number of edges and contains one cycle only. Taking 
into account that each element Xj is represented by a cyclic graph with 
2q' + 1 edges, the only possibility to get rid of the "extra" cycles in 5k and is 
to glue the elements Xj all long one of each other. To get maximal number 
of edges in 5k, we keep the cyclic structure of Xj. When doing this, we 
have to choose the direction and the edge of Ai to glue the first elements of 
Aj, j = 2, . . . , A;. The number of possibilities is given by 2*^-1 X (2g' + 1)*^-^ 
and we get the first equality of (6.13). 

Regarding the case of even q = 2q' , we conclude that to construct a tree 
with the maximal number of edges we have first to construct a tree from 
each element Xj and then to draw A; — 1 arcs between these k trees obtained. 
It is easy to see that the maximal size of the tree obtained by gluing the 
edges of Xj is q' and the number of different trees is given by Catalan number 
{q' + ) [25]. Then we produce the connected diagrams exactly as in 
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the case of dilute Y model. This gives the second relation of (6.13). Lemma 
is proved. 



Similarly to the F-model, the sparse random graphs regime of X-model 
is more complicated than the dilute random graphs regime. The average 
value of random variable Xn'^^ = Tryl^^' is studied in [1, 13] (see also 

[16]). The limit ^(3) = lim„^oo ^^^^n determined by a system of 

recurrent relations that have the form similar to that described in Lemma 

6.2. We do not present these results. 

(2a') 

The variance of Xk^ ' is studied in [26] in more general setting than that 
of this paper. The resulting expression is also determined by a system of 
recurrent relations. This system is much more complicated that that for the 
first moment of We refer the reader to the paper [26] for corresponding 
theorems. 



6.5 Formal limit of the partition function 

Returning to the normalized partition function (2.12) that describes the 
y-model with q = 2, one can see that 

PnU^ 2p„n \ 1 + e Pn-n^ >- J 



(6.15) 

where (3' = f) — tgn and g^ = {pnn) ^. It is clear that the three asymp- 
totic regimes introduced in Section 4 are characterized by the corresponding 
behavior of the inverse "temperature" /? and elementary analysis shows that 



n - 1 , / 1 + ~ 

hm 



' 0, if /? = Const; 

0, if /3 = ilogf , K c< n; 

exp{2tc-i}-i ^ if /? = i log ^, c = Const, 



where pn = e ^'^(l + e ^^) ^. Then, assuming that the limit of the last 
term of (6.15) exists, we conclude that (4.4) implies relation 

hm -l^logZr.iP,tgn) = %^ (exp{^ 1 - l) + £ ^^^T^' (6-16) 



where u indicates the full, the dilute, and the sparse random graphs regimes, 
and Soj,3 denotes the Kronecker (5-symbol. 

Comparing expressions (6.1), (6.2), and (6.4) for with uj = 1,2,3, 
we see the difference between the limiting rate functions of (6.16). Indeed, 
does not depend on c in dilute random graphs regime, oj = 2, and is 
given by a polynomial in degrees of 1/c in the sparse random graph regime, 
a; = 3. In particular, regarding the small-t expansion of the right hand-side 
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of (6.16) in the sparse random graphs limit and taking into account that 
(3) = 1 + c ^, we see that this approximation to the rate function is 
t(l + f)(l + o(l)),t^O. The same difference between the rate functions is 
observed in the general case of with exponential term of (6.16) replaced 
by i(exp{^} - 1). 

Let us stress that the above reasoning relies strongly on the existence of 
the limit 

lim -\ log E;3' {e*^"^" } = fS) (6.18) 

with Vn = Xn^ or Yn'^\ This problem is difficult to solve by using the 
cumulant expansion (3.1). As it is pointed out many times (see for exam- 
ple [8]), the terms of the formal relation (3.1) with positive g are used to 
enumerate combinatorial structures related to the matrix integrals. In the 
rigorous sense, the series (3.1) never converges in this case. Existence and 
analyticity in g of the terms of asymptotic expansion (3.1) with respect to 
n is proved for the matrix integrals just recently [3, 7] by using powerful 
techniques of the integrable models and the Riemann-Hilbert problem. 



7 Summary 

The Gibbs weight determined by the graph Laplacian generates a measure 
fin on the ensemble of n-dimensional adjacency matrices of simple non- 
oriented graphs. This measure is invariant with respect to the permutations 
of the basis vectors and determines the Erdos-Renyi ensemble of random 
graphs with the edge probability p„. 

Regarding the sum over the set of weighted adjacency matrices as the 
analog of the matrix integrals, we determine the discrete analog of matrix 
models with the quartic potential. In the general case of g-power potential, 
we distinguish two different families of discrete Erdos-Renyi matrix models. 
These are related with the numbers of g-stcp walks and g-step closed walks 
over the random graphs denoted by Yn'^^ and Xn \ respectively. 

The logarithm of the standard partition function of these models is de- 
termined by the cumulants of random variables Xn^ and Y^^^ . We develop a 
diagram technique to study the limiting behavior of the cumulants in three 
major asymptotic regimes of full, dilute and sparse random graphs deter- 
mined by the edge probability as n — >■ oo. We prove that the limits of 
these cumulants, when properly normalized, exist in all of the three asymp- 
totic regimes. As a consequence, the Central Limit Theorem is shown to be 
true for centered and normalized variables X and Y. This implies CLT for 
the moments of the normalize spectral measure of the adjacency matrix of 
random Erdos-Renyi graphs. 

We show that the limiting expressions of the cumulants are related with 
the number of non-rooted trees constructed with the help of labeled edges. 
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In the simplest case of the dilute random graphs regime, the exponential 
generating function of these numbers H{z) verifies the Lagrange (or Polya) 
equation. Passing to the sparse random graphs regime, we derive more 
general equations that determine H{z) in this case. 

These results imply an observation that the asymptotic regimes we con- 
sider are different not only with respect to the normalization factors of the 
cumulants, but also with the respect to the rate functions of large deviations 
formally determined as the limit fq{t) (6.18). It should be noted that we did 
not prove rigorously the existence of this limit because our results concern 
the leading terms of the cumulants only. 
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